DO THE QCD SUM RULES SUPPORT FOUR-QUARK STATES? 
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We test the validity of the QCD sum rules applied to the light scalar mesons, the charmed mesons 
D s o(2317) and D s i(2460), and the X(3872) axial meson, considered as tetraquark states. We find 
that, with the studied currents, it is possible to find an acceptable Borel window only for the 
X(3872) meson. In such a Borel window we have simultaneouly a good OPE convergence and a 
pole contribution which is bigger than the continuum contribution. We interpret these results as 
a strong argument against the assignment of a tetraquark structure for the light scalars and the 
D s0 (2317) and D s i(2460) mesons. 
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I. INTRODUCTION 

From the light scalar mesons to the heavy 
"chamonium-like" A (3872), there are now many 
states that do not fit comfortably in the spectrum of 
constituent quark model predictions. The light scalar 
states below 1.5 GeV are too numerous to be accommo- 
dated in a single qq multiplet and the nature of these 
states has been a source of controversy for over 30 years 
[![. The lightest nonet is composed, in principle, by the 
isoscalars cr(600) and /o(980), the isodoublet k(800) and 
the isovector ao(980). In a naive qq assignment it is hard 
to explain the /o - ao mass degeneracy and why a and 
k are broader than the other two. The strange-charmed 
mesons £> s0 (2317) and D s i(2460) (J p = 0+and 1+ 
respectively) are too light to fit in the quark model 
prediction, with the Dm(2317) lying about 160 MeV 
below most predictions 0. The AT (3872), with quantum 
numbers J PC = 1 ++ , does not fit in the charmonium 
spectrum and presents a strong isospin violating decay, 
disfavoring a cc assignment [2j . 

The structure of all these states has been exten- 
sively discussed and many alternatives have been pro- 
posed: meson molecules, four-quark states, glueballs 
(in the case of scalars) and hybrids {qgq). The idea 
that the light scalar mesons could be four-quark bound 
states has been first proposed by Jaffe in 1977 Q, and 
has later been extrapolated to heavier sectors. Jaffe 
proposed that some states may be composed of two 
quarks and two antiquarks (qqqq) arranged so that the 
(anti)quark-(anti)quark correlation is important, forming 
what is called a (anti)diquark. Recently the existence of 
tctraquarks received some support from lattice calcula- 
tions [4|, which, however, are not yet definitive. 

The QCD sum rules (QCDSR) 0,0,0 have been pre- 
viously used to study the light scalars [3, 0, [T3, E| , the 
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strange-charmed scalars [H, EH and the A(3872) [3| as 
diquark-antidiquark states. 

In Q it was assumed that the light scalars were 
tetraquarks and no attempt to compute their masses in 
QCDSR was performed. Instead a calculation of their 
decay widths, using their experimental masses was pre- 
sented. At the same time, in [l2[ the masses of several 
charmed scalars were calculated. With the tetraquark 
hypothesis, the decay width of the D s j (2317) and of the 
A (3872) were calculated in [ll| and in [16[ respectively. 

While the masses were often very close to the experi- 
mental values, the widths were not always as narrow as 
found in experiments. This is expected because, unless 
some symmetry violation is involved, tetraquarks can de- 
cay more easily since no quark pair creation is needed and 
a "fall - apart" decay is allowed. 

From 2003 to 2006, the QCDSR calculations of masses 
and decay widths evolved rapidly and became much more 
rigorous. While the first calculations aimed only at es- 
timating some order of magnitude and only the Borel 
stability was checked, the last ones were much more con- 
cerned with OPE convergence and with pole dominance, 
which are traditional tests, from which one can determine 
the quality of the calculation. 

The improvement of the standards was also motivated 
by "the pentaquark experience" . In this case from the 
begining there was an experimental controversy about 
the very existence of the particle. After the first round 
of promising results, it was realized [17U18| that the pen- 
taquark sum rules were problematic, because it was al- 
ways very difficult to find a Borel window in which one 
would have at the same time good OPE convergence and 
pole dominance. In favor of the QCDSR practitioners it 
must be said that sum rules with more than three quarks 
presents new and challeging aspects. The number of 
possible interpolating currents increases significantly and 
also one has to worry about subtracting the two-hadron- 
reducible contributions [T^ | , a problem never encountered 
before in QCDSR calculations. Finally, to make things 
even more complex, there may be a mixing between two 
and four-quark states. This requires the combination of 
interpolating fields of different dimensions with the in- 



2 



traduction of a new parameter. 

Relating pentaquarks and tetraquarks may be very in- 
structive. In both cases negative results were gradu- 
ally found, but there was always still a lot of work to 
be done, such as computing higher order contributions 
to the OPE, instanton contributions, a s corrections and 
new possible interpolating currents. Therefore it took a 
long time until a negative opinion about pentaquarks was 
formed in the QCDSR community. We have now gath- 
ered evidence to believe that QCD sum rules calculations 
of tetraquark properties have reached the same turning 
point found before in the case of pentaquarks. This is 
the point where, even though there are improvements to 
be made, we do not believe that these improvements will 
change the conclusion of a series of works pointing to the 
non-existence of tetraquarks. 

In this work we review some of the tetraquark sum 
rules with special attention to the validity limits of the 
method. In section II we work out the sum rules of 
the axial strange-charmed Z? s i(2460) as a prototype for 
this analysis, and extend the discussion to other states. 
The study of the D s i(2460) complements the calculations 
published in [l2j • In section III we extend the analysis of 
section I to the light scalars, studying some of the inter- 
polating fields proposed for these states and study also 
the charmed scalars. In section IV we examine the sum 
rules for the X(3872). 



II. THE CHARMED AXIAL MESON D s i(2460) 

The interpolating operator for D s i(2460) (as a 
diquark-antidiquark state) is built by extension of the 
operator used to describe D s0 (2317) in ref. [12|], chang- 
ing the diquarks so we get an axial current: 



V2 



[(ulC lb c b ){u dllJl Cs T e ) + u^d], (1) 



where a, 6, c, ... are color indices and C is the charge 
conjugation matrix. We choose to work with an axial 
light antidiquark to avoid instanton contributions to the 
sum rule [2(| • 

The sum rule for the charmed axial meson is con- 
structed from the two-point correlation function: 

H, v {q) = i [d*x e^(0|T[i M (x)jt(0)]|0> = 



= -Ri{<i 2 )(g f ^<i 2 - iniv) + ^o(q 2 )q^- (2) 

Since the axial vector current is not conserved, the two 
functions, IT and IIo, appearing in Eq. ^fy are indepen- 
dent and have respectively the quantum numbers of the 
spin 1 and mesons. 

The calculation of the phenomenological side proceeds 
by inserting intermediate states for the axial vector me- 
son and parametrizing its coupling to the current j^, in 
Eq. {l}, m terms of the meson decay constant fo sl as: 



<0|j M |£> s i) = y/2f. 



A>i TO D sl £ M 



(3) 



the phenomenological side of Eq. ([2]) can be written as 

up;r{q2) = ijk^ t + w\ + . . . f (4) 



where the Lorentz structure <? M „ projects out the spin 1 
state. The dots denote higher axial- vector resonance con- 
tributions that will be parametrized, as usual, through 
the introduction of a continuum threshold parameter sq 

In the OPE side we work at leading order and con- 
sider condensates up to dimension six. We deal with the 
strange quark as a light one and consider the diagrams 
up to order m s . To keep the charm quark mass finite, 
we use the momentum-space expression for the charm 
quark propagator. We calculate the light quark part of 
the correlation function in the coordinate-space, which is 
then Fourier transformed to the momentum space in D 
dimensions. The resulting light-quark part is combined 
with the charm-quark part before it is dimensionally reg- 
ularized at D = 4. 

We can write the g^ structure of the correlation func- 
tion in the OPE side in terms of a dispersion relation: 

P(s) 



ds- 



r 



(5) 



where the spectral density is given by the imaginary part 
of the correlation function: p(s) — ^Im[H OPE (s)]. After 
making a Borel transform on both sides, and transferring 
the continuum contribution to the OPE side, the sum rule 
for the g^v structure can be written as 

- 2fl 3i ml si e~™ 2 sl* 12 = / ds e-^ 12 p(s) , (6) 



where p(s) = p pert (s) + p^(s) + p {G ' 2) (s) + p mlx {s) + 
p { 1 q)2 (s), with 
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(3 + a), (7) 
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6m s 4(qq) 
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(*) = 
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p mlx {s) 



2- OH ' j ^ 

m s (&{qgo.Gq) + (sga.Gs)(2 - 3n ) j + 
2m ^- Gq \l-oa + 2a 



(10) 
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< 5g >= (s) = M| fdaK^s), (11) 
12tt^ J a 

where A = m 2 c / s and K{a, s) — m 2 — as. For the charm 
quark propagator with two gluons attached we used the 
momentum-space expressions given in ref. Q . 

In order to extract the mass m£, sl without knowing 
about the value of the decay constant /d s1 , we take the 
derivative of Eq. ^ with respect to 1/M 2 , divide the 
result by Eq. ([6]) and obtain: 

2 Q ds e-V^ s pjs) 
m °«- Qdse-^pis) ■ {U) 

In the numerical analysis of the sum rules, the val- 
ues used for the quark masses and condensates are [22j : 
m c = 1.23 GeV, m s = 0.1 GeV, (qq) = -(0.23) 3 GeV i , 
(qga.Gq) = m 2 (qq) with [7] m 2 = 0.8 GeV 2 and 
(g 2 G 2 ) = 0.88 GeV 4 . We evaluate the sum rules for 
three values of sq: y/so — 2.7 GeV, ^/sq — 2.9 GeV and 
3.1 GeV. 



A. Pole versus continuum 

We get an upper limit for M 2 by imposing that the 
QCD continuum contribution should be smaller than the 
pole contribution. The maximum value of M 2 for which 
this constraint is satisfied depends on the value of sq. The 
comparison between pole and continuum contributions 
for y/s~o = 2.9 GeV is shown in Fig. [TJ 

100. s 1 ' 2 = 2.9GeV 



■ - Pole 
Continuum 
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FIG. 1: The dashed line shows the relative pole contribution 
(the pole contribution divided by the total, pole plus contin- 
uum, contribution) and the solid line shows the relative con- 
tinuum contribution. The pole contribution should be bigger 
than the continuum, which happens for M 2 < 1.5 GeV 2 for 
^/sE = 2.9 GeV. 



The same analysis for the other values of the contin- 
uum threshold gives M 2 < 1.4 GeV 2 for v /s7^ = 2.7 GeV 
and M 2 < 1.6 GeV 2 for y/sE = 3.1 GeV. 

In Fig. [21 we show the D s i mass obtained from 
Eq. (JT2J) , in the M 2 region below the upper limit obtained 
above. We limit ourselves to the region M 2 > 1.2 GeV 2 
where the curves are more stable. Averaging the mass 
over all this region we get: 

2.7-, 
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FIG. 2: The D„i mass as a function of the sum rule parameter 
(M 2 ) for different values of the continuum threshold. The 
arrows indicate the region allowed by the upper limit imposed 
by the dominance of the QCD pole contribution. 

m Dal = (2.3 ±0.2) GeV , (13) 

which is compatible with the experimental value 
£ s i(2460) ||. 

B. OPE convergence 

There is however a stronger constraint to the lower 
bound of the M 2 region. We have to analyze the conver- 
gence of the OPE by comparing the relative contribution 
of each term in Eqs. ([7]) to (fTTj) . to the right hand side 
of Eq. (O. The series converges better for higher values 
of M 2 , so that requiring a good convergence sets a lower 
limit to M 2 . This analysis in shown in figure [3l 

Figure[3]shows no convergence in any region allowed by 
the upper bound given by pole/continuun analysis. This 
means that the lower bound given by OPE convergence 
will be higher than the upper bound, and there is no 
"sum rule window" where we can completely trust the 
results for this current. 

The results above illustrate very well how we can re- 
produce the mass of a given state and then after a more 
carefult analysis conclude that the state is not a parti- 
cle as such, being rather one of the possible continuum 
excitations. 
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FIG. 3: The OPE convergence in the region 1.0 < M < 
1.7 GeV 2 for y / 5o' = 2.9 GeV. We start with the perturbative 
contribution and each subsequent line represents the addition 
of a condensate of higher dimension in the expansion. 



III. THE SCALAR MESONS 

A. Light scalars 

The same situation described in the last section is en- 
countered in many sum rules with interpolating operators 
built with more than three quark fields. The light scalar 
meson interpolating operators used in ref. [9J] are: 

Cafcc (-dec Cj 5 d b ) (ildj5 C¥ e ) , 



3a 


— ^abc^dec 




^abc^dec 


3fo 


V2 


3a 


tabc^dec 


V2 


3k 


— ^abc^dec 



e a bc(dec{UaCj 5 d b )(q d j5Cs^), q = u,d. (14) 

They yield very low upper limits to M 2 when the pole and 
continuum contributions are analysed: M 2 < 0.73 GeV 2 
for a (980) and / (980), M 2 < 0.62 GeV 2 for the k(800) 
and M 2 < 0.54 GeV 2 for the cr(600). The analysis was 
performed with the same parameters used in 0: Sq = 
1.0 GeV 2 , = 1.2 GeV 2 , 4° = 1-5 GeV 2 . 

In figure [4] we show the OPE convergence for ao (980) 
and /o(980) (which have the same sum rule), in the same 
way shown in figure [3] In this figure we see that there 
is no OPE convergence in any region allowed by the up- 
per bound. In fact, the situation of the light scalars is 
even worse than that of the Z3 s i(2460), since the rela- 
tive contribution of the dimension-6 condensate is even 
bigger. A possible reason for this is the fact that we 
are working with very small values for the Borel Mass 
(M 2 < 1 GeV 2 ). As a matter of fact, once the integral 
on the right hand side of Eq.[5]is evaluated, the OPE side 
becomes a series with decreasing powers of Af 2 , which 



FIG. 4: The a (980)//o(980) OPE convergence in the region 
0.4 <M 2 < 0.7 GeV 2 for y/sE = 1-2 GeV. 



eventually become negative so that higher condensates 
will be divided by higher and higher powers of M 2 . In 
the case of the tetraquarks the series begins with M 10 and 
one still has a positive power of M 2 for the dimension-8 
condensate. However, it is hardly justifiable to truncate 
the series at this point since higher dimension conden- 
sates will be proportional to (1/M 2 ) 2 f where D is the 
dimension of the condensate and for M 2 < 1 GeV 2 these 
condensates will not be suppressed, at least for D ~ 10. 

It is interesting to notice that the authors of ref. [23| 
have arrived at the conclusion that the ao(980) scalar 
meson is not a four-quark state using a different criterion. 
The authors of ref. [23J have annalyzed the QCD sum 
rules of the do (980) meson considered as a normal two- 
quark state, and also as a four-quark state. While they 
could reproduced both the mass and width of the ao (980) 
considered as a two-quark state, they were not able to 
reproduce the width of the ao(980) considered as a four- 
quark state. 

It could be argued that these problems are related with 
the specific currents that we are working with, and that 
there could be other currents that might work better. In 
Ref. [IH, five different interpolating operators for each 
of the light scalar mesons have been tested. In the case 
of the a these currents were: 



SI = 

V° = 

n = 

A a - 

P? = 



(ulCj 5 d b )(u a j 5 CdJ - u b j 5 CdJ), 

(ulc^Muar^cdJ - u b ^ 75 CdZ), 

{u T a CG^d h ){u a a» v CdJ - u h a» v CdF a ), 

[u T a c^d b ){u a ^cdJ - u bl »cdZ), 

{u T a Cd b ){u a CdJ -u b CdZ). (15) 



The currents for the other light scalars can be obtained 
by the following substitutions: k : (ud)(ud) — > (ud)(ds), 
fo : (ud)(ud) — > (us)(us) + (ds)(ds) and ao : [ud)(ud) — > 
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(us)(us) — (ds)(ds). The authors of [Tl| have tested 
all these currents and various linear combinations and 
found out that the better results were obtained with the 
particular combination: r/^ — cos9Aq + smBV^, with 
cos0 = l/\/2. They also obtain good results for the other 
light scalars with similar combinations. 

We used the same analysis used above with the spectral 
densities obtained in ll| and agree that the OPE con- 
vergence up to dimension 8 is quite good. On the other 
hand the pole dominance requirement imposes very low 
upper limits to M 2 : M 2 < 0.8 GeV 2 for a or f (^/sE = 



IV. HEAVIER TETRAQUARKS 

The situation improves as the quarks in the interpolat- 
ing operator become heavier. In the case of the V(3872), 
the following operator was used in ref. [l4[ : 



■ X ^^-abc^dec 

3 p. 



{llC-f b c b ){q d -i^Ccl) 



V2 

(llCj^iqd^CcJ)}. 



(17) 



The continuum contribution analysis for sets the up 
1.6 GeV), M 2 < 0.45 GeV 2 for k (y^o" = 1.2GeV)and per limit at M 2 < 2.6 GeV 2 for a threshold of 
M 2 < 0.35 GeV 2 for a (yio = 1. GeV), 

This means that the whole sum rule window lies below 
M 2 < 1 GeV 2 and, as commented above, it is at least 
dangerous to truncate the series at this order. 



4.3 GeV. The OPE convergence in this region is shown 
in figure [S We see that, for M 2 > 1.9 GeV 2 , the addi- 



B. Charmed scalars 

In the case of the charmed scalar D s0 (2317), the cur- 
rent used for it in ref 11211 is: 



tabc^dec 



^/2 



[{u T a C lb c b ){u dlb Cs T e ) + u <-> d] . (16) 



If we require that the pole contribution be bigger than 
the continuum contribution we obtain M 2 < 1.37 GeV 2 
for ^/so = 2.7 GeV. In figure [5] we show the OPE con- 
vergence for the current (fl"6|) and we see that the OPE is 
still not convergent in the allowed region, as in the case 
of tme meson D sl (2460) studied in the previous section. 
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FIG. 5: The j a OPE convergence in the region 1.0 < M < 
2.0 GeV 2 for ^ = 2.7 GeV. 
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FIG. 6: The j* OPE convergence in the region 1.5 < M 2 < 
3.1 GeV 2 for = 4.3 GeV. 



tion of a subsequent term of the expansion brings the 
curve (representing the sum) closer to an asymptotic 
value (which was normalized to 1). Furthermore the 
changes in this curve become smaller with increasing di- 
mension. These are the requirements for convergence 
and in this case we get a sum rule window in the re- 
gion 1.9 GeV 2 < M 2 < 2.6 GeV 2 . The mass obtained in 
[14| considering the allowed Borel window is 



m x = (3.92 ±0.13) GeV, 



(18) 



which is compatible with the experimental value 
V(3872). 

A similar situation is found if we replace the c quarks 
in Eq. (jTTJ) by b quarks in order to predict the X b mass 
(as done in ref. [HI)- ^ n this case the allowed Borel 
window is in the region 6.0 GeV 2 < M 2 < 7.0 GeV 2 , and 
the predicted mass is 

m Xb = (10.14 ± 0.11) GeV , (19) 

which is in agreement with the findings in ref. (24| . 
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From what was seen above we can conclude that for 
heavier tetraquarks the sum rules satisfy the validity cri- 
teria and hence allow the determination of the masses 
of these states. However, even in the present case we 
can not yet be very positive. Firstly because, as usual, 
the calculations might still be improved, with, for exam- 
ple, the inclusion of a s corrections. Secondly because it 
remains very difficult to reproduce the X narrow decay 
width, as shown in [16]. If the A(3872) is proved to be 
a tetraquark state, it still remains to explain why we do 
not observe tetraquark states with charge different from 
cc states, such as (cu)(cd) or (cd)(cu) states, which would 
also have trustable QCDSR as the A (3872). In this sense, 
the observation of a double charmed meson ((cc)(qq)), 
which sum rule also obey all the convergence and pole 
dominance criteria [25], would be very important to re- 
ally determine the existence of tetraquark states. 

V. CONCLUSION 

We have performed a QCD sum rules calculation of 
the -D s i(2460) mass considering this state as a tetraquark 
and reanalized other recent similar tetraquark sum rules, 
giving special attention to the validity criteria of the 
method. We found that in the case of the lighter states, 
a(600), k(800), a (980), / (980) and also in the case of 
the intermediate -D s n(2317) and D s i(2460) states, for the 
currents used in refs. 0, [l2| . there are no values of the 



parameters sq and M 2 that satisfy all the desired condi- 
tions. In order to obtain results from the sum rules for 
these states we must abandon one or more of the condi- 
tions and choose the parameters arbitrarily. 

When the interpolating operator is constructed with 
heavier quark fields the situation becomes better. We 
found suitable regions for the A (3872) and its extension 
to the bottonic sector the Xf,. 

This problem was also present in the case of the pen- 
taquarks [l7| and seems connected to the high dimension 
of many-quark states interpolating operators, indepen- 
dently of the exact form of these operators. This may be 
an indication from the sum rules that light many-quark 
states can not be considered as ressonances separated 
from the continuum. Heavier many-quark states are sup- 
ported by the sum rules in what concerns their masses. 
However it is very difficult (if possible) to explain their 
narrow decay widths. 

While one might always argue that the so far existing 
calculations could be improved and the final conclusions 
might still change, to us at this point in time, this seems 
unlikely. 



Acknowledgements 

This work has been partly supported by FAPESP and 
CNPq. 



[1] C. Amsler and N. Tornqvist, Phys. Rept. 389, 61 (2004), 
and references therein. 

[2] E. S. Swanson, Phys. Rept. 429, 243 (2006) and refer- [11 
ences therein. [12 

[3] R. L. Jaffe, Ph ys. Rev. D15, 267 (1977); R. L. Jaffe, 

|hep-ph/000Tl23| [13 

[4] H. Suganuma, K. Tsumura, N. Ishii and F. Okiharu, [14 
PoS LAT2005, 070 (2006); |hep-lat/0509121 ; F. Oki- 
haru, H. Suganuma and T. T. Takahashi, Phys. Rev. D [15 
72, 014505 (2005). [16 

[5] M.A. Shifman, A.I. and Vainshtein and V.I. Zakharov, 

Nucl. Phys. B147, 385 (1979). [17 

[6] L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rep. 

127, 1 (1985). [18 

[7] For a review and references to original works, see e.g., 

S. Narison, QCD as a theory of hadrons, Cambridge [19 
Monogr. Part. P hys. Nucl. Phys. Cosmol. 17, 1 (2002) 
hep-h/0205006 ; QCD spectral sum rules , World Sci. 
Led. Notes Phys. 26, 1 (1989); Acta Phys. Pol. 26 [20 
(1995) 687; Riv. Nuov. Cim. 10N2 (1987) 1; Phys. 
Rept. 84, 263 (1982). [21 

[8] J. I. Latorre and P. Pascual, Jour. Phys. Gil, L231 

(1985); S. Narison, Phys. Lett. B175, 88 (1986). [22 

[9] T. V. Brito, F. S. Navarra, M. Nielsen, M. E. Bracco, [23 
Phys. Lett. B608, 69 (2005). [24 
[10] H-J. Lee and N.I. Kochelev, Phys. Lett. B642, 358 

(2006); Z.-G. Wang, W.-M. Yang, S.-L. Wan, J [25] 
Phys. G31, 971 (2005); A. Zhang, Phys. Rev. D61, 



Zhang, T. Huang and T. Steele, 



114021 (2000); A. 
|hep-ph/0612146"l 
H.X. Chen, A. Hosaka and S.L. Zhu, |hep-p h/0609163 
M. E. Bracco, A. Lozea, R. D. Matheus, F. S. Navarra, 
M. Nielsen, Phys. Lett. B624, 217 (2005). 
H. Kim and Y. Oh, Phys. Rev. D72, 074012 (2005). 
R. D. Matheus, S. Narison, M. Nielsen, J. M. Richard, 
Phys. Rev. D75, 014005 (2007). 
M. Nielsen, Phys. Lett. B 634, 35 (2006). 
F. S. Navarra and M. Nielsen, Phys. Lett. B 639, 272 
(2006). 

R. Matheus and S. Narison, Nucl. Phys (Proc. Suppl.) 
152, 263 (2006). 

R.D. Matheus, F. S. Navarra and M. Nielsen, Braz. J. 
Phys. 36, 1397 (2006). 

Y. Kondo, O. Morimatsu and T. Nishikawa, Phys. Lett. 
B 611, 93 (2005); Y. Kwon, A. Hosaka and S. H. Lee, 
|hep-ph/050504u] 

H-J. Lee and N.I. Kochelev, Phys. Lett. B642, 358 
(2006). 

B.L. Ioffe, Nucl. Phys. B188, 317 (1981); B191, 591(E) 
(1981). 

S. Narison, Phys. Lett. B466, 345 (1999). 
CA. Domingues and N. Paver, Z. Phys. C39, 39 (1988). 
F.-K. Guo, P.-N. Shen, H.-C. Chiang and R.-G. Ping, 
Nucl. Phys. A761, 269 (2005). 

F.S. Navarra, M. Nielsen and SH. Lee, |hep-ph/0703071| 



